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We report a rich spectrum of isolated solitons residing inside [embedded into) the continuous radi- 
ation spectrum in a simple model of three- wave spatial interaction in a second-harmonic-generating 
planar optical waveguide equipped with a quasi-one-dimensional Bragg grating. An infinite sequence 
of fundamental embedded solitons are found, which differ by the number of internal oscillations. 
Branches of these zero-walkoff spatial solitons give rise, through bifurcations, to several secondary 
branches of walking solitons. The structure of the bifurcating branches suggests a multistable con- 
figuration of spatial optical solitons, which may find straightforward applications for all-optical 
switching. 



I. INTRODUCTION 

Recent studies have revealed a novel class of embedded 
solitons (ESs) in various nonlinear-wave systems. An ES 
is a solitary wave which exists despite having its inter- 
nal frequency in resonance with linear (radiation) waves. 
ESs may exist as codimension-one solutions, i.e., at dis- 
crete values of the frequency, provided that the spectrum 
of the corresponding linearized system has (at least) two 
branches, one corresponding to exponentially localized 
solutions, the other one to delocalized radiation modes. 
In such systems, quasilocalized solutions (or "generalized 
solitary waves" ) in the form of a solitary wave resting 
on top of a small-amplitude continuous-wave (cw) back- 
ground are generic ||]. However, at some special values of 
the internal frequency, the amplitude of the background 
may exactly vanish, giving rise to an isolated soliton em- 
bedded into the continuous spectrum. 

Examples of ESs are available in water-wave mod- 
els, taking into account capillarity and in several 
nonlinear-optical systems, including a Bragg grating 
incorporating wave-propagation terms Q and second- 
harmonic generation in the presence of the self-defocusing 
Kerr nonlinearity (the latter model with competing 
nonlinearities was introduced earlier in a different con- 
text [§). 

It is relevant to stress that ESs, although they are iso- 
lated solutions, are not structurally unstable. Indeed, 
a small change of the model's parameters will slightly 
change the location of ES (e.g., its energy and momen- 
tum, see below), but will not destroy it, which is clearly 
demonstrated by the already published results [^|J|]. In 
this respect, they may be called generic solutions of codi- 
mension one. 

ESs are interesting because they naturally appear 



when higher-order (singular) perturbations are added 
to the system, which may completely change its soli- 
ton spectrum. Optical ESs have a potential for applica- 
tions, due to the very fact that they are isolated solitons, 
rather than occurring in continuous families. The stabil- 
ity problem for ESs was solved in a fairly general analyti- 
cal form in Rcf. Q , which was also verified by direct sim- 
ulations of the model considered. It was demonstrated 
that ES is a semi-stable object which is fully stable to 
linear approximation, but is subject to a slowly growing 
(sub-exponential) one-sided nonlinear instability. Devel- 
opment of this weak instability depends on values of the 
system's parameters; in some cases, it is developing so 
slowly that ES, to all practical purposes, may be regarded 
as a fully stable object 0. 

In the previously studied models, only a few branches 
of ESs were found, and only after careful numerical 
searching, which suggest they may be hard to observe in 
a real experiment. The present work aims to investigate 
ESs in a recently introduced model of a three-wave inter- 
action in a quadratically nonlinear planar waveguide with 
a quasi-one-dimensional Bragg grating |^ , which can be 
quite easily fabricated. It will be found that ESs occur in 
abundance in this model, hence it may be much easier to 
observe them experimentally. It should also be stressed 
that, unlike the previously studied models, in which ESs 
appear in relatively exotic conditions, e.g., as a result of 
adding singular perturbations |^ or specially combining 
different nonlinearities the model that will be con- 
sidered below and found to give rise to a rich variety of 
ESs, is exactly the same which was known to support 
vast families of ordinary (non-embedded) gap solitons. 
This, in particular, implies that ES can be found in the 
corresponding system under the same conditions which 
are necessary for the observation of the regular solitons, 
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i.e., the experiment may be quite straightforward. An es- 
timate of the relevant physical parameters will be given 
at the end of the paper. 

The rest of the paper is organized as follows. In sec- 
tion 2, we recapitulate the model and obtain solutions in 
the form of fundamental zero-walkojf ^Ss, which, physi- 
cally, correspond to the case when the Poynting vector of 
the carrier waves is aligned with the propagation direc- 
tion. The analysis is extended in section 3 to the case of 
fundamental walking ESs, for which the Poynting vector 
and the propagation distance are disaligned. Concluding 
remarks are collected in section 4. 

II. THE MODEL AND ZERO-WALKOFF 
SOLITONS 

The model describes spatial solitons produced by the 
second-harmonic generation (SHG) in a planar waveg- 
uide, in which two components of the fundamental har- 
monic (FH) , vi and V2 , are linearly coupled by the Bragg 
reflection on a grating in the form of a system of scores 
parallel to the propagation direction z (for a more de- 
tailed description of the model see [|| ) : 

K'Vl)z+i{vi)x+V2+V'iV2=0, (1) 
i{v2)z-i{v2)x+Vi+V'iVl=0, (2) 

'^i{v'i)z-qv3. + D{v:>,)xx+viV2=Q. (3) 

Here is the second-harmonic (SH) field, a; is a normal- 
ized transverse coordinate, q is a real phase-mismatch 
parameter, and D is an effective diffraction coefficient. 
The diffraction terms in the FH equations ([l]) and (||) 
are neglected as they are much weaker than the artificial 
diffraction induced by the Bragg scattering, while the SH 
wave, propagating parallel to the grating, undergoes no 
reflection, hence the diffraction term is kept in Eq. (||). 

Experimental techniques for generation and observa- 
tion of spatial solitons in planar waveguides are now well 
elaborated ( |^), and the waveguide carrying a set of 
parallel scores with a spacing commensurate to the light 
wavelength (which is necessary to realize the resonant 
Bragg scattering) can be easily fabricated. Therefore, the 
present system provides a medium in which experimen- 
tal observation of ESs is most plausible. As mentioned 
above, the observation of ES in this system should be fur- 
ther facilitated by the fact that it supports a multitude 
of distinct ES states, see below. 

Eqs. (l)-(3) have three dynamical invariants: the 
Hamiltonian, which will not be used below, the energy 
flux (norm) 

/ + 00 
[\vi{x)\^ + \V2{X)\'' + ^v:i\^]dx , (4) 
-oo 

and the momentum. 



/ + 00 
{{v,)lv, + {v2)lv2 + 2{v,)lv:i)dx. (5) 
-oo 

The norm played a crucial role in the analysis of the ES 
stability carried out in 1^. 




FIG. 1. The {k,q) parameter plane of the three- wave 
model (0) - (|^. The linear analysis (the results of which 
are summarized in the inset boxes) shows that ES with c = 
may occur only in the region between the solid bold lines. The 
bundle of curves emanating from the point (k = l,g = —4) 
are branches of embedded-soliton solutions with c = 0. 

Soliton solutions to Eqs. (l)-(3) are sought in the form 

vi,2{x, z) = exp(iA:z) mi,2(?), ^^3(2^, z) = exp(2ifc2;) , 

(6) 

where ^ = x — cz, with c being the walkojf (slope) of the 
spatial soliton's axis relative to the light propagation di- 
rection z. The substitution of (|^) into Eqs. (1)~(3) leads 
to an 8th-order system of ordinary differential equations 
(ODEs) for the real and imaginary parts of wi.2,3 (primes 
standing for d/dS): 

- kui + i{l - c)u[ + U2 + u-iU2 — 0, (7) 

~kU2 — j(I + c)?i2 + Ml + ?i3M^ — 0, (8) 

— (4fc + q)u^ + -DU3 — 2«cu3 + U1U2 = 0. (9) 

Before looking for ES solutions to the full nonlinear 
equations, it is necessary to investigate the eigenvalues A 
of their linearized version, in order to isolate the region 
in which ESs may exist. Substituting ui,U2 ~ exp(A^), 
and it3 ~ exp(2A^) into Eqs. (6)~(8) and linearizing, 
one flnds that the FH and SH equations decouple in the 
linearized approximation. The FH equations give rise to 
a biquadratic characteristic equation, 

(1 - c^fX" + 2 [(1 + c2)yfc2 _ (1 _ c^)] A2 -f (fc2 -1)2 = 0, 

(10) 
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and the SH equation produces another four eigenvalues 
given by 



[DA^ - {4k + q)Y + 4c'^X^ =0. 



(11) 



A necessary condition for the existence of ESs is that 
the eigenvalues given by Eq. ( p^ ) have non-zero real parts 
- this is necessary for the existence of exponentially lo- 
calized solutions - while the eigenvalues from Eq. ( p^ ) 
should be purely imaginary (otherwise, one will have or- 
dinary, rather than embedded, solitons). This discrimi- 
nation between the two sets of the eigenvalues is due to 
the fact that Eqs. (0) and (^) for the FH components 
are always linearizable, while the SH equation may 
be nonlinearizable, which opens the possibility for the ex- 
istence of ESs 1^ . As it follows from Eqs. dl^) and du]), 
these two conditions imply 



k^ + c^ < 1: 



Ak + q< c^/D. 



(12) 



For the case c = 0, the parametric region defined by the 
inequalities ( [l^ ) is displayed in Fig. 1. 



certain double-humped two-solitons (bound states of two 
fundamental solitons) appear to be stable [||). Note that 
double- and multi-humped ESs must exist too as per a 
theorem from Ref. [^, but leaving them aside, we will 
still find a rich structure of fundamental ESs. 




FIG. 3. A diagram of the c = embedded solitons on the 
(energy-flux, mismatch) plane. The inset zooms the most in- 
teresting part of the diagram. 




FIG. 2. Typical examples of the fundamental embedded 
solitons with the zero walkofF: (a) the ground-state for fc = 0; 
(b) the same solution for k — —0.95; (c,d) the first and eighth 
"excited states" for k — 0. 



In Ref. numerous ordinary {gap |10|| ) soliton solu- 
tions to the present model have been found by means of a 
numerical shooting method. To construct ES solutions, 
we applied the same method to Eqs. (0), (|) and (|), 
allowing just one parameter to vary. From each ES so- 
lution that was found this way, branches of the solutions 
were continued in the parameters k, q and c, by means of 
the software package AUTO [|ll| . Note that the c = so- 
lutions admit an invariant reduction U2 = — u*, M3 — u^, 
which reduces the system to a 4th-order ODE system, 
thus making numerical shooting feasible. 

We confine the analysis to fundamental solitons, which 
implies that the SH component U3 has a single-humped 
shape (a distinctive feature of gap solitons in the same 
system is that not only fundamental solitons, but also 



We begin with a description of the results from the re- 
duced case c = 0, when an additional scaling allows us 
to set D = 1 without loss of generality. The results are 
displayed in Figs. 1-3. There is a strong evidence for ex- 
istence of an infinite "fan" of fundamental ES branches. 
Among them, we define a ground-state soliton as the one 
which has the simplest internal structure (Fig. 2a). The 
next "first excited state" differs by adding one (spatial) 
oscillation to the FH field (Fig. 2c). Adding each time 
an extra oscillation, we obtain an indefinitely large num- 
ber of "excited states" (as an example, see the 8th state 
in Fig. 2d). We stress, however, that all the "excited 
states" belong to the class of the fundamental solitons, 
rather than being bound states thereof. 

In Fig. 1, the first nine states (branches) are shown in 
the (fc, q) parametric plane. Note that the whole bundle 
of the branches originates from the point {k = l,q = —4), 
which is precisely the intersection of the two lines which 
limit the existence region of ES (see Eq. ( p^ with c — 0). 
At this degenerate point, the linearization (see above) 
gives four zero eigenvalues. More branches than those 
depicted in Fig. 1 have been found, the numerical re- 
sults clearly pointing towards the existence of infinitely 
many branches, accumulating on the border q + Ak — 
of the ES-region. In the accumulation process, each U3 
component is successively wider, while the ui^2 ones have 
more and more internal oscillations. 

Since k is an arbitrary propagation constant, from 
physical grounds, the results obtained for the c = solu- 
tions are better summarized in terms of energy flux E vs. 
mismatch q (Fig. 3). Note that all the branches shown in 
Fig. 3 really terminate at their edge points, which exactly 
correspond to hitting the boundary fc = — 1, see Fig. 1. 



3 



It is also noteworthy that all the solutions are exponen- 
tially localized, except at the edge point fc = — 1, where 
a straightforward consideration of Eqs. (0)-(||) demon- 
strates that, in this case, ES are weakly (algebraically) 
localized as \x\ ^ oo (cf. Fig 2b): 



ui « ^-{4k + q)\x\-\, U2 «(l/2)v/-(4fc + 9)1x1^2^ 

Finally, we observe from Figs. 1 and 3 that the first 
"excited-state" branch has a remarkable property that it 
corresponds to a nearly constant value of q. This means 
that while, generally, ES are isolated [codimension-one) 
solutions for fixed values of the physical parameters, this 
branch is nearly generic, existing in a narrow interval of 
the q-values between —4.0 and —3.74. 



III. WALKING SOLITONS 

We now turn to ESs with c 7^ 0, i.e., walking ones. 
These were sought for systematically by returning to the 
full 8th-order-0DE model and allowing the AUTO pack- 
age to detect bifurcations (of the pitchfork type), while 
moving along branches of the c = solutions. It tran- 
spires that all the bifurcating branches have c ^ 0, i.e., 
they are walking ESs. Such solutions are of codiniension- 
two in the parameter space (i.e., the solutions can be 
represented by curves k(q), c{q)), which can be estab- 
lished by a simple counting argument after noting that 
the 8th-order linear system has two pairs of pure imagi- 
nary eigenvalues. Alternatively, the walking ESs can be 
represented, in terms of the energy flux and momentum 
(see Eqs. (^) and (||)), by curves E{q) and P{q). We 
present results only for the walking solutions which bi- 
furcate from the ground and first excited c = states, 
while other walking ESs can also be readily found. 

It was found that the ground-state branch has ex- 
actly two bifurcation points, giving rise to two distinct 
walking-ES solution branches (up to a symmetry) . These 
new branches are shown, in terms of the most physi- 
cally representative c{q) and E{q) dependences, in Fig. 
4. Note that they, eventually, coalesce and disappear. As 
the inset to Fig. 4b shows, they disappear via a tangent 
(fold, or saddle-node) bifurcation. 

The first excited state has three bifurcation points. 
One of them gives rise to a short branch of walking ESs 
that terminates, while two others appear to extend to 
q = —00 (their ostensible "merger" in Fig. 5 is an arti- 
fact of plotting). It is known that, in the large- mismatch 
limit q ^ —00, the present three- wave model with the 
quadratic nonlinearity goes over into a modified Thirring 
model with cubic nonlinear terms [ p^ . This suggests that 
the latter model may also support ES. However, consid- 
eration of this issue is beyond the scope of the present 
work. 





(a) 








FIG. 4. Two branches of "walking" (c 7^ 0) embedded soli- 
tons bifurcating from the ground-state c — branch: (a) the 
walkoff c, and (b) the energy flux E vs. the mismatch q. The 
horizontal segment in (a) shows the branch of the c = solu- 
tions. The inset in (b) shows that the two branches meet and 
disappear via a typical tangent bifurcation. 

Fig. 4 clearly shows that, in a certain interval of the 
mismatch parameter q, the system gives rise to a multi- 
stability, i.e., coexistence of different types of spatial soli- 
tons in the planar optical waveguide (for instance, taking 
account of the fact that each c ^ branch has symmetric 
parts with the opposite values of c, we conclude that there 
are five coexisting solutions at q taking values between 
about —8 and —11). This situation is of obvious interest 
for applications, especially in terms of all-optical switch- 
ing Indeed, switching from a state with a larger value 
of the energy flux to a neighboring one with a smaller flux 
can be easily initiated by a small localized perturbation, 
in view of the above-mentioned one-sided semistability 
of ES, shown in a general form in Such a switching 
perturbation can be readily made controllable and mov- 
able if created by a laser beam launched normally to the 
planar waveguide and focused at a necessary spot on its 
surface ||l^. Switching between the two branches with 
c 7^ can be quite easy to realized too, due to small 
energy-flux and walkoff/momentum differences between 
them, see Fig. 4. 



IV. CONCLUSION 

To conclude the analysis, it is necessary to estimate the 
actual size of the relevant physical parameters. This is, 
in fact, quite easy to do, as there is no essential difference 
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in the estimate from that which was presented in Ref. ||] 
for the ordinary sohtons in exactly the same model. This 
means that a diffraction length ^ 1 cm is expected for 
the SH component, and, definitely, the diffraction lengths 
for the FH components, which are subject to the strong 
Bragg scattering, will be no larger than that. Thus, a 
sample with a size of a few cm may be sufficient for the 
experimental observation of ESs. The sample may be an 
ordinary planar quadratically nonlinear waveguide with a 
set of parallel scores written on it. The other parameters, 
such as the power of the laser beam that generates the 
solitons, etc., are expected to be the same as in the usual 
experiments with the spatial solitons As concerns 
the weak semi-instability of ESs, it may be of no practi- 
cal consequence for the experiment, as it would manifest 
itself only in a much larger sample. In this connection, it 
may be relevant to mention that, strictly speaking, the 
usual spatial solitons observed in numerous experiments 
are all unstable (e.g., against transverse perturbations) in 
the usual (linear) sense, but the instability has no room 
to develop in real experimental samples. 

Finally, we see from Figs. 4 and 5 that the maximum 
walkoff that ESs can achieve is, in the present notation, 
slightly smaller than 1. According to the geometric inter- 
pretation of the underlying equations (||) - (||) (see details 
in the original work Q), this implies that the maximum 
size of the misalignment angle between the propagation 
direction and the axis of the spatial soliton may be nearly 
the same as the (small) angle between the Poynting vec- 
tors of the two FH waves and that of the SH wave. 

To summarize the work, we have found a rich spec- 
trum of isolated solitons residing inside the continuous 
spectrum in a simple model of the three- wave spatial in- 
teraction in a second-harmonic-generating planar optical 
waveguide equipped with a quasi-one-dimensional Bragg 
grating. An infinite sequence of fundamental embedded 
solitons were found. They differ by the number of in- 
ternal oscillations. The embedded solitons are localized 
exponentially, except for a limiting degenerate case, when 
they become algebraically localized. Branches of the 
zero-walkoff spatial solitons give rise, through bifurca- 
tions, to several branches of walking solitons. The struc- 
ture of the bifurcating branches provides for a multistable 
configuration of the spatial optical solitons. This may 
find straightforward applications to all-optical switching. 
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FIG. 5. Three branches of the walking (c 7^ 0) embedded 
solitons bifurcating from the c = branch corresponding to 
the first "excited state" , depicted similarly to Fig. 4 The inset 
in (a) shows in detail the central part of the diagram. 
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